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^ ■ Abstract 
G^ 

Oh| We compute the ratio A^/Ajj^ between the scale parameter A/,, associated 

■ with a lattice formulation of QCD using the overlap-Dirac operator, and Aj^^ 

■ of the MS renormalization scheme. To this end, the necessary one-loop rela- 
CN , tion between the lattice coupling qq and the coupling renormalized in the MS 

scheme is calculated, using the lattice background field technique. 
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I. INTRODUCTION 



The overlap-Dirac operator that has been derived from the overlap formulation of 
chiral fermions on the lattice provides a lattice regularization of massless QCD, i.e. of 
^ ■ chiral fermions coupled vectorially to a gauge field, without the need of fine tuning. It 
preserves chiral symmetry without fermion doubling, circumventing the Nielsen-Ninomiya 
^ ! theorem |^. The overlap-Dirac operator Do of a massless fermion can be written as (a is 



H ' the lattice spacing) 



Do = ^[l + j,0], (1.1) 



where O depends on the link variables, and is Hermitian with eigenvalues ±1. The simplest 
such example is given by the Neuberger-Dirac operator [|l| 



= - c 
a 



1 + X{X^X)-^/^], (1.2) 



X = Dw - -P, (1.3) 

a 

where is the Wilson-Dirac operator (with the Wilson parameter r set to its standard 
value, r = 1) 
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1 



w 
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(V; + V J - aV^V J , (1.4) 
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V^ip{x) = -[U {x, fi)ip{x + afi) — ip{x)] , (1.5) 

and p is a real parameter subject to the constraint < p < 2.Q The normahzation constant 
c = p can be absorbed in the definition of the fermion fields. Unlike the Wilson-Dirac 
operator D^^i, is not analytic in the link variables when the operator X in Eq. ( |1.3| ) has 
a zero eigenvalue. As discussed in Ref. 0], this is the price to pay for putting strictly massless 
fermions on the lattice. However, such a lack of analiticity is expected to be harmless in the 
continuum limit P,|23[|. 

Djss satisfies the Ginsparg- Wilson relation ^ 

^,D + D^, = aD^,D, (1.6) 



which protects the quark masses from additive renormalization [|1],0. The Ginsparg- Wilson 
relation allows us to write, at finite lattice spacing, relations that are essentially equivalent 
to those holding in the low-energy phenomenology associated with chiral symmetry (see e.g. 
Refs. P,^). It indeed implies the existence of an exact chiral symmetry of the lattice action 
under the transformation |T0[ 



5ij{x) = 75(1 - |aD)^(a;), <5V'(a;) = ^(a;)(l - |aD)75, (l-T) 

thus leading to chiral Ward identities that ensure the non-renormalization of vector and fia- 
vor non-singlet axial vector currents, and the absence of mixing among operators in different 
chiral representations. The axial anomaly then arises from the non-invariance of the fermion 
integral measure under flavour- singlet chiral transformations [|T^-[TB|. It is also worth men- 
tioning that lattice gauge theories with Ginsparg- Wilson fermions have been proved to be 
renormalizable to all orders of perturbation theory |T^]. 



The important point is that lattice Dirac operators satisfying Eq. (|1.6|) are not affected 
by the Nielsen-Ninomiya theorem P], thus they need not suffer from fermion doubling. A 
lattice formulation of QCD satisfying the Ginsparg- Wilson relation would overcome the 
complications of the standard approach (e.g. Wilson fermions), where chiral symmetry is 
violated at the scale of the lattice spacing. As a consequence of chiral symmetry, leading 
scaling corrections are O(a^), as opposed to 0{a) of the chiral symmetry violating case. 

Indeed, D-^ avoids fermion doubling. However, its locality properties in the presence 
of a gauge fleld are not obvious. is not strictly local, and locality should be recov- 
ered only in a more general sense, i.e. allowing an exponential decay of the kernel of 
with a rate which scales with the lattice spacing and not with the physical quantities. In 
Ref. P3[ the locality of Z)n has been proved for sufficiently smooth gauge flelds. Moreover 



numerical evidence has been presented for typical gauge flelds in present-day simulations. 
Thus, the Neuberger-Dirac operator seems to have all the right properties that a lattice 
Dirac operator should have in order to describe massless quarks. A major open question 
seems to be its practical implementation, since appears much more complicated than 



Nonperturbatively one expects —mc < p < 2, where rric is the negative critical mass associated 
with the Wilson-Dirac operator. 
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the usual Wilson-like Dirac operators. In this respect some progress has been achieved (see, 
e.g., Refs. |17,23,18-p0|,^,pl-25[|), and simulations may become feasible in the near future. 

In this paper we calculate the ratio between the A-parameter associated with a lattice 
formulation of QCD using the overlap-Dirac operator (A^) and that of the MS renormal- 
ization scheme (Aj^). In order to have a complete discretization of QCD, we will consider 
the Wilson formulation for the pure gauge part of the theory. Actually, the part of the cal- 
culation involving fermions is independent of the regularization chosen for the gluonic part. 
We recall that the ratio of any renormalization group invariant quantity to the appropriate 
power of Ai approaches a constant in the continuum limit go ^ 0. Indeed A^ is a particular 
solution of the renormalization group equation 



d d ^ 



Al = 0(a) 



I.e. 



a A] 



exp 



90 dg 



f^Lig) 
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1 + O {gt) 
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where PL{go) is the lattice /5-function, and bo,bi the first two coefficients of its perturbative 
expansion: 



bo 
bi 



(47r)2 V 3 3 ^ 



(47r) 



34 



-N 



3 N 



:i.io) 

:i.ii) 



in SU{N) gauge theory with Nf fermion species. The calculation of the ratio Al/Ajj^ 
requires a one-loop perturbative calculation on the lattice. As it turns out, the necessary 
computation is much more cumbersome than in the case of Wilson fermions, due to the 
more complicate structure of the Neuberger-Dirac operator. 



II. FORMULATION OF THE PROBLEM 

The lattice regularization of QCD we consider is described by the action 

= - E Tr [1 - u,,{x)] + M^)D^{^, y)4^^iy)■ (2.1) 

9o x,fi,u i=l x,y 

where U^^{x) is the usual product of link variables U^{x) along the perimeter of a plaquette 
originating at x in the positive fi-u directions, and Nf is the number of massless fermions 
considered. We recall that A^^ is independent of the fermionic masses, so the results of our 
calculation will also hold for the massive cases. f\ 



In order to describe massive fermions one may write the overlap-Dirac operator in the form 1 26,19] 
Do oc 1 -|- ;U -|- (1 — ^)^^0. Values < < 1 describe fermions with mass m-j, < mf < oo. For 
small m/ is proportional to /U. 
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In order to evaluate the ratio Al/Aj^ we need to calculate, at one loop order, the relation 
between the lattice coupling go and the renormalized coupling g of the MS scheme: 



go = Zg{go,aij)g, 

where /i indicates a renormalization scale. Writing 

Zg{go, x)2 = 1 + g^ (26o lnx + lo) + 0{g*o), 

one has 

lo = 26oln (Al/Ams)- 



(2.2) 



(2.3) 



(2.4) 



The computation of Zg is easier in the background field gauge. In fact, this renormaliza- 
tion constant has a simple relationship with the background field renormalization constant 



Za 27 



ZA{go,afi)Zg{go,afif 



1. 



(2.5) 



As a consequence of this relation, in order to calculate Zg one only needs to calculate the 
one-loop self-energy of the background field. 

The MS renormalized one-particle irreducible two-point function of the background field 
is given by 



AA/^\ab 



rr(p) 



\ i=i 



'^'n'ip/i^) /g' 



(2.6) 



where 



N 
16^ 



11, p2 205 3 1 



36 2A 4A2 



+ 



2 

3^V 
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(2.7) 



and A is the gauge parameter. On the lattice one writes 



i-E^?oV«(«p) 



i=l 



hi 



(2.8) 



The bare and renormalized functions are related by 



i-E^?oV«(«p) 



i=l 



(2.9) 



Therefore, using relation (|2.5|), 



i-E.=i^gM^nM 

' l-E,,=i^7^V«(p//x)' 



(2.10) 



So, in order to evaluate /q we have to perform the one-loop calculation of the function 
z/'^^)(ap) in the background field gauge as formulated on the lattice P5|,p5[. We mention that 
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the background field technique on the lattice has been recently employed for the calcula- 
tions of the third coefficient of the lattice /5-function for Wilson-like lattice formulations of 
QCD [pO|-p^. We refer to the above cited references for the relevant formulae. 

In order to perform the lattice perturbative calculation of u^^\ap) we must formally 
expand in powers of go- The weak coupling expansion of has been already discussed 
in Ref. |3^. We list here the relevant formulae for our one- loop calculation. 

Let us first write down the weak coupling expansion of the Wilson-Dirac operator D^^. 
This will be useful for constructing the relevant vertices of • We write 



= Dw(g,p) - -P = Xoip)i2n)''5\q-p)+X,{q,p)+X2iq,p) + 0{g-^), (2.11) 
a 



where 



i 1 1 
= - ^li, sin ap^ + - ^{l- cos ap^) p, 



(2.12) 



Xi{q,p) =90 J d'k6{q -p- k)A,{k)Vi,^{p + k/2), 
ViA^) = i'lf.cosaq^ + sinag^, 



(2.13) 



X2{q,p) 



Qq f d^ki d'^k2 



6{q-p-ki- k2)A^{ki)A^{k2)V2AP + ki/2 + k2/2), (2.14) 



2 J (27r)4 
^2,M(g) = -i^f^asinaq^ + a cos aq^. 

Let us write the Fourier transform of the Neuberger-Dirac operator in the form 
D^{q,p) = D,{p){2TxY5\q - p) + S(g,p). 
Dq{p) is the tree level inverse propagator: 

° 2Kp) + 6(p)] +2' 



where 



UJ 



(P) = ~ ( ^ ^^^^ + E*^^ ~ cosap^) - pY 



1/2 



Kp) = - 11(1 - cosap^) - -p. 



One can easily check that for < p < 2 the propagator has only one pole, at p 



(2.15) 



(2.16) 



(2.17) 
(2.18) 
0. The 



function 'E{q,p) can be expanded in powers of go as |33 
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1 2 

FIG. 1. Diagrams with fermionic lines contributing to the one-loop function i'^^\ap). Dashed 
lines ending on a cross represent background gluons. Solid lines represent fermions. 



aE{q,p) 



1 
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Xo{p)Xl{p,q)Xo{q) 



X 



(27r)4 uj{p) + uj{q) u;{p) + uj{k) u{q) + uj{k) 
-Xoip)Xl{p, k)X,{k, q) - X,{p, k)Xl{k)X^{k, q) - X,{p, k)Xl{k, q)X^{q) 
uj{p) + u){q) + uo{k) 



-Xoip)Xlip, k)Xoik)Xl{k, g)Xo(g) 



+ .. 



2.19) 



uj{p)uj{q)uj{k) 

From one can read off the vertices necessary for the one loop calculation of u^^^ap). 



III. RESULTS AND DISCUSSION 

Two diagrams containing fermions contribute to u^^^p), shown in Figure 1. Given that 
the 4-point vertex contains a part with an internal momentum {k in Eq. ( p.l9| )), the corre- 
sponding part of the second diagram actually has the same connectivity as the first diagram. 

The algebra involving lattice quantities was performed using a symbolic manipulation 
package which we have developed in Mathematica. For the purposes of the present work, 
this package was augmented to include the propagator and vertices of the overlap action. 

The one-loop amplitude u^^^ap) can be written as: 



z/(i)(ap) = iy^^>{ap)\Mf=o + NfJ2 

i 

(the index i runs over the two diagrams shown in Figure 1), where, 

'' In a^p^ ^ 



(3.1) 



op h'i{ap) = do,i + a p < di^i + d2 



(47r) 



0{{apY 



(3.2) 
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TABLE I. Coefficients 



p 




di,2 


0.2 


-0.00222139 


0.01803836 


0.3 


-0.00418268 


0.01753315 


0.4 


-0.00563107 


0.01732216 


0.5 


-0.00681650 


0.01727870 


0.6 


-0.00785477 


0.01736058 


0.7 


-0.00881237 


0.01755678 


0.8 


-0.00973458 


0.01787211 


0.9 


-0.01065738 


0.01832255 


1.0 


-0.01161402 


0.01893460 


1.1 


-0.01263969 


0.01974719 


1 9 


— u.UlOl / uou 


n 09081 fifin 


1.3 


-0.01507869 


0.02222125 


1.4 


-0.01662489 


0.02408181 


1.5 


-0.01853492 


0.02658162 


1.6 


-0.02101054 


0.03002396 


1.7 


-0.02443278 


0.03495886 


1.8 


-0.02966205 


0.04255372 



(p2 _ 4 sin^(p^/2)). The coefficients dn^i depend on p, but not on N or Nj . Note that 
the diagrams in Figure 1 do not involve the gluonic propagator, so they are independent of 
the choice of regularization for the pure gluonic part of the action. 

To extract the p-dependence, we ffist isolate the divergent terms; these are responsible 
for the logarithms. There are only a few such terms, and in the pure gluonic case their 
values are well known. We can use these values also in diagrams with fermions, applying 
successive subtractions of the type: 




(3.3) 



where is the inverse fermionic propagator. All remaining terms now contain no diver- 
gences, and can be evaluated by Taylor expansion in ap. 

At this stage, one is left with expressions which no longer contain p and must be nu- 
merically integrated over the loop momentum. Given the complicated form of the overlap 
vertices, these expressions turn out to be quite lengthy, containing a few thousand terms. 

The integration is done in momentum space over finite lattices; an extrapolation to 



infinite size is then performed, in the manner of Ref . . We evaluated the integrals for a 
range of values of the parameter p, as presented in Table 1. For all values of p tabulated, 
lattice sizes L < 128 were sufficient to yield answers with at least 5 significant digits (the 
uncertainty coming from a systematic error in the extrapolation, which can be estimated 
quite accurately). As the endpoints of the perturbative domain of p are approached (p — *• 
0, p — > 2), increasingly larger lattices are required for similar accuracy; this is, of course, 
a reflection of the divergences in the propagator at these endpoints. As shown in Figure 2, 
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FIG. 2. A^/A^ in ^[/(S), as a function of the parameter p. Nf = (dotted line), A''^ 
(dashed hne), A'^^ = 3 (sohd hne). 



the value of the A^-parameter varies also in a more pronounced manner near the endpoints, 
where it is expected to diverge. 

The coefficients dn,i obey two constraints, which we have used as verifications of our 
procedure: 

i) Gauge invariance requires J2i do,i = 0. One may check this property both on the 
algebraic expressions and on the numerical results. Substituting the numerical results 
for the coefficients, we find zero within the error estimates, for all values of p ; this 
serves also to verify our estimation of systematic errors. 

ii) The coefficients d2^i must coincide with those of the continuum. We checked that this 
is so: (i2,i = 2/3 , ^2,2 = . 



Thus, our result takes the form 

'^^^\ap) |ao=i= i^^^\ap) |jv^=o, Ao=i +Nf 



247r2 



In(aV) + kfip) 



(3.4) 



where 
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z/(^)(ap) Itv =0, Ao=i = 



lliV_ 2 2^ 1 



In(aV) - ^ + 0.21709849Ar 



48^2 - ' 



(3.5) 



is the pure gauge part ||3^ , |35| , |30| , and = (ii i + rfi 2 is presented in Table 1. We note 

again that the choice of the pure gluonic part of the action affects the first, but not the 
second, summand in Eq. (|3.4|) . 

From Eqs. (O, O, pJOT), we find 



/o = 0.16995599Ar + A^. 

8A^ 



727r^ 



The values of A^/Aj^^g follow immediately from Eq. (|2] 
interest are (N = 3): 



^L/-^Ms\Nf=0 
^L/^Ms\Nf=l, p=l 
^L/-^Ms\Nf=3, p=l 



0.034711 
0.025042 
0.011273 



(3.6) 

Some particular cases of 

(3.7) 
(3.8) 
(3.9) 



For comparison we report the corresponding values of the ratio A^/A^ for the standard 
Wilson formulation of fermions 13611 



^L/^Ms\Nf=l, r=l 
^I-/ ^Ms\Nf=3, r=l 



0.029412 
0.019618. 



(3.10) 
(3.11) 
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